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We study local diffusion rate D in Chirikov standard map near the critical golden curve. Numerical 
simulations confirm the predicted exponent q = 5 for the power law decay of D as approaching 
the golden curve via principal resonances with period qn {D ~ l/q^)- The universal self-similar 
structure of diffusion between principal resonances is demonstrated and it is shown that resonances 
of other type play also an important role. 
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During last years intensive investigations have allowed 
to understand the structure of critical invariant curves 
at chaos border in Hamiltonian dynamical systems with 
divided phase space |^,^. Usually the analysis is carried 
out for two dimensional (2d) area preserving maps, the 
paradigm being Chirikov standard map |^. The critical 
invariant curve is characterized by the rotation number 
and its continued fraction expansion. It has been shown 
that the small scale structure near invariant curves with 
the same tail in this expansion is universal for all smooth 
Hamiltonians with two degrees of freedom and for 2d 
maps Among all invariant curves better studied 

are those with golden rotation number rg = (v^ — l)/2, 
whose expansion is a series of I’s. This rg is the most 
irrational number and therefore it is believed that invari¬ 
ant curves with golden tails are locally the most robust 
ones. The structure of the critical golden curve has been 
studied by means of renormalization group approach and 
it has been shown that the phase space structure is self¬ 
similar and universal on small scales. 

Different scaling exponents have been found in this 
critical regime and they were successfully used to de¬ 
termine the diffusion rate through the destroyed invari¬ 
ant curve called cantorus However, no flux passes 

through the golden curve at the critical value of per¬ 
turbation parameter K = Kg. In this case the tra¬ 
jectory has only local diffusion rate D, which depends 
on its distance Ar„ to the golden curve. This diffu¬ 
sion characterizes the motion of a particle in the vicinity 
of invariant curve Vg at different levels n of convergents 
Tn = Pn/On of the Continued fraction expansion of Vg 
(n = I/I,r 2 = 2/3, rs = 3/5...). One can expect a 
power law dependence of D on the resonant approximant 
g„, namely D ~ qn~°‘■ Chirikov gave a simple argument 
for a = 5 1^. According to him D ~ (Arji)^/tjj, where 
Ar„ = |rg — r„| ~ and is the typical inverse 

frequency of small oscillations around principal reso¬ 
nance qn- Then, On, ~ qnAujn, where Aujn is the width 
of resonance qn 1^. In the critical case Chirikov overlap 
criterion implies Atu„ ^ Ar„, which gives 


D « ADo/ql ~ (Ar„)5/2 ^ 


( 1 ) 


where Dq = is the quasilinear diffusion rate 

and 5yn = \yn — yg\ is the distance of unstable periodic 
orbit 2/„ with rotation number to the golden curve yg 
along the symmetry line. The exponent v = 2.14699... 
can be determined from the exponent a for Syn l/^Zn 
found in H] (j^ = 5/cr). 

The fast decay of diffusion rate near the chaos border 
rg means that a diffusing particle will never reach the bor¬ 
der itself. This slow diffusion gives a long sticking of tra¬ 
jectory around stable islands on different renormalization 
levels. As a result the statistics of Poincare recurrences P 
(integrated probability to return into a given region after 
a time larger than r) decays with r as P{t) ~ I/t^. Such 
kind of decay has been first observed in [Q , where the av¬ 
erage value of the exponent p ~ 1.5 has been found. Fur¬ 
ther investigations have shown that the power law decay 
of P is a generic property of Hamiltonian systems with 
divided phase space [Hj|]. However, according to numeri¬ 
cal results the power itself is not universal, varying in the 
range 1 < p < 2. Moreover different maps with golden 
chaos boundary give different p values , which seems to 
be in contradiction with universal self-similar structure 
of phase space near the golden border. Indeed renor¬ 
malization arguments give p = 3 0, which is in sharp 
contradiction with numerical results. One of the possible 
reasons of the above contradiction is sticking of particles 
near the stable islands between principal resonances qn- 
However, the attempts to take into account these inter¬ 
mediate resonances gave p = 2 | W|, which is too large 
compared to the numerical results for the golden border 
1 < p < 1.35 [^. Therefore the problem of Poincare 
recurrences remains unsolved and more detailed inves¬ 
tigations of the phase space structure near the golden 
border should be performed. One of the reasons why the 
properties of P{t) are so important is that the correlation 
function of dynamical variables C (r) and the probability 
p{t) to stay in a given region for a time t > t are related 
to P{t) 
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C'(r) ~ /i(r) ~ tP{t)/ < t >~ 1/t^ ^ (2) 

where < t > is the average return time. The above rela¬ 
tions follow from the ergodicity of motion on the chaotic 
component of the phase space |Q,|[P. The decay of corre¬ 
lations with power pc = p— 1 < 1 can lead to a divergence 
of global diffusion rate and to strong fluctuations in di¬ 
vided phase space. 

In this Letter we investigate the behavior of local diffu¬ 
sion rate D near the critical golden curve. Our first aim 
was to verify the theoretical prediction and to analyze 
the structure of D at different renormalization levels. As 
a model we have chosen Chirikov standard map 

y = y — iL/(27r) sin(27ra;) , x = x + y modi (3) 

with perturbation parameter corresponding to critical 
golden curve K = Kg = .97163540631.... To mea¬ 
sure the local diffusion rate D = (Ay)^/Af we apply 
the efficient method used for the investigations of Arnold 
and modulational diffusion in Ref. {t is number of 
iterations). This method allows to measure very small 
diffusion rates (down to computer noise level) with rel¬ 
atively small number of iterations. By this method, 
we compute D at different resonant Fibonacci approx- 
imants r„ = Pn/qn of the critical curve Vg. We use 
from Np = 10 to Np = 100 trajectories near unstable 
periodic orbits of period y„ (these points had been de¬ 
termined by MacKay |^). Each trajectory is integrated 
for about T = 1000 x iterations. The total interval 
T is divided into Ny^ = 10 windows , where the aver¬ 
age y displacement was computed with the smoothing 
function / = (ntNyi/T). Usually we take /3 = 4;6. 
Such smoothing allows to suppress regular oscillations by 
a factor oc (Au,/T)^^+^ [^. To control the accuracy of 
our numerical computation of D, we determine D also 
near stable periodic orbits in the center of resonance (/„, 
where our method gives a value D ~ 10“^^, which corre¬ 
sponds to the level of computer round-off error in double 
precision. 

Our results in Fig. 1 confirm the theoretical predic¬ 
tion (U) for the variation of D over more than 20 orders of 
magnitude. The numerical fit gives a = 4.99 ± 0.02 and 
A = 0.0066. We attribute the fluctuations at small 
values to the fact that for a large number of iterations tra¬ 
jectories can exit from the chaotic layer corresponding to 
the initial q„. This effect disappears for larger where 
the local diffusion rate is sufficiently small or for shorter 
number of iterations. Let us note that the numerical 
value of A is surprisingly small. Our explanation for this 
fact is the following. According to the action change 
after half period of rotation in the chaotic separatrix layer 
is quite small Ay ^ 4A^ exp(—7rA/2) (here A = 27r is the 
frequencies ratio). This gives an order of magnitude esti¬ 
mate for A = D{qn = l)/73o ~ (^2/)^/(2Do) ~ 0.003 and 
explains its small value. A more accurate estimate of A 
requires taking into account higher orders of perturbation 


in K. The measured diffusion rate in the chaotic com¬ 
ponent is well separated from the diffusion in the stable 
regions produced by numerical round-off errors (Fig.l). 

The above result shows the global structure of the dif¬ 
fusion rate D while approaching the chaos border via 
resonant approximants qn- However, an interesting ques¬ 
tion concerns the behavior of D in between = Pn/qn 
and Tn+i ■ The comparison of D on these scales should re¬ 
flect the self-similar structure of phase-space on different 
renormalization levels. To check this self-similarity we 
measure D on two symmetry lines x = 0; 0.5. The sym¬ 
metry line a; = 0.5 crosses the main part of chaotic layers 
and contains mainly unstable points, while the other line 
X = 0 passes mainly through stable islands. The known 
structure of periodic orbits on symmetry lines of map (H) 
1^ allows to find the basic renormalization intervals on 
stable X = 0 and unstable x = 0.5 lines. 

For the unstable line the first renormalization level in¬ 
terval is Ayi = |yi — yyj where y„ is the y value of the 
unstable periodic orbit with rotation number on the 
line X = 0.5. The next interval Ay 2 = |y 4 — yio| lies on 
the other side of the golden curve. The m—th interval is 
= |y 3 m -2 —y 3 m-i- 4 |- This m—series selects the subse- 
quence of n—values which we will denote by Um = 3m—2. 
Intervals with odd values of m lie above the invariant 
curve Tg and those with even m lie below. The self¬ 
similarity of the phase space implies that the dependence 
of diffusion D on the position inside m—th interval should 
be approximately the same as in to -|- 1—th interval after 
rescaling. At larger to— values the self-similarity is 
expected to become better and better. To check numer¬ 
ically this self-similarity we compute the diffusion rate 
at 320 homogeneously distributed points in the intervals 
Aym for TO = 1,..., 6. To compare different intervals we 
rescale the diffusion rate defining Dr = where 

Dm is the diffusion at to— th level. The position inside 
each interval is denoted by Ay/j = |(y —y 3 m- 2 )|/Aym for 
y between y 3 m -2 and y 3 m-i -4 (0 < Ay^; < 1). 

The comparison of renormalized diffusion rate Dr on 
two levels to = 4 and to = 6 on the unstable line is 
shown in Fig. 2. The diffusion rate is self-similar in 
agreement with universal phase space structure near the 
golden curve. The minimal diffusion rate Dr on lev¬ 
els TO = 4; 6 is determined by computer noise and is 
different on the two levels due to the different normal¬ 
ization factors. The rare fluctuations in the upper dif¬ 
fusion plateau are presumably due to exits of trajec¬ 
tory from initial chaotic layer. Points with high diffu¬ 
sion rate are in the chaotic component, while those with 
minimal diffusion correspond to trajectories in stable is¬ 
lands. The self-similarity was also observed at other to 
levels, both when intervals were on the same side or on 
opposite sides of the golden curve. This self-similarity 
becomes better with the growth of to. The sharp sep¬ 
aration between two levels of diffusion allows to deter¬ 
mine the leading resonances in each renormalization in- 
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terval. The biggest gap in the diffusion is for peri¬ 
odic orbit (labelled b in Fig. 2) with rotation number 
rt = 1,2,1], where the triplet of I’s in curly 

brackets is repeated to — 1 times. This rotation num¬ 
ber is not from the series of principal resonances given 
by the continued fraction expansion of Vg and labelled 
by i in Fig. 2 (n = 1,1,1,1,1]). While reso¬ 

nance rb is significantly larger than the principal one 
there are also other resonances which are of comparable 
or smaller size than r^: rc = 1, 2,1,1,1,1], 

rd = [{lll}^_i,l,2,2,2,l], re = 1, 2, 2,1], 

Tg = [{lll}^_i,l,l,l,3,2,l], rb = 1,1,1,- 

2,1,1,1], n = [{lll}^_i,l,l,l,l,l,2,l,l]. It is inter¬ 
esting to remark that most of these rotation numbers 
have continued fraction expansions containing mainly I’s 
and 2’s, which is in qualitative agreement with a conjec¬ 
ture made in on the basis of analysis of the Schmidt- 
Bialek fractal diagram of map (||). This fact is also in 
agreement with numerical results in ||l^ . 

In Fig. 2 we see also other types of resonances, namely 
ra = Pa/qa = [{11 l}^_i, 1,4,1] and rj = pf/qf = 
[{111}^_^, 1,1,1, 3, Ij. However, they have a different 
structure than previous ones, displaying two and three 
intersections with unstable line x = 0.5 giving rise to 
double (a) and triple (/) drops in diffusion. Indeed the 
periodic orbit inside resonance (r/) has period 2qa 
(3(7/). These orbits are not present in the limit K ^ 0, 
and correspond to a new chain of islands in the chaotic 
layer around resonances Va {r/)- The important conse¬ 
quence of the analysis of the structure of these resonances 
is that, in spite of the fact that renormalization group de¬ 
scribes quite well the convergence to golden curve, we see 
that on each renormalization level resonances other than 
those of the main series of Vg, and even some which are 
not present at IF = 0, occupy a sizeable part of phase 
space. The existence of such non-standard resonances 
might explain the lack of universality of the exponent p in 
(^ and the numerical value of p significantly less than 2. 
Indeed, very long Poincare recurrences can be originated 
by sticking of orbits not only near the main resonances 
but also around these non-standard resonances and the 
chains of islands around them. The general description 
of the phase-space structure should take into account the 
presence of these resonances. An interesting question is 
which is the weakest cantorus on a given renormalization 
level which will determine the transition time between 
different renormalization levels. According to Fig.2 all 
drops in D are associated to periodic orbits and not to 
invariant curves. This means that at our accuracy level 
there is no other invariant curve except r^. 

We have also studied the local diffusion on the sta¬ 
ble line X = 0.0 (Fig. 3). The size of renormalization 
interval n is defined as Aj/„ = \yn — yn+ 2 \ and the vari¬ 
able Ayn = \{y — yn)\/Ayn. The structure of diffusion 
is also self-similar. As expected, the main part of the 
renormalization interval is occupied by stable islands. 


The largest resonances with rotation number pa = [{!}„] 
and Pi = [{1}^,1,1] correspond to the main series of 
Tg. There are also other resonances: pc = [{l}jj,3,1], 
Pd = [{!}„, 3,1,1], Pe = [{!}„, 2,1], P/ = [{!}„, 2,1,1,1], 
Pg = 1,1], Ph = [{!}„, 2, 2,1]. Again the biggest 

resonances have only I’s and 2’s in the continued frac¬ 
tion. Some of the resonances seen on a; = 0.0 are also 
observed on a; = 0.5 and one can easily establish the 
correspondence. The resonance pb has the same rotation 
number as pa, but it corresponds to a different orbit. The 
period of this orbit is 9 times larger than that of pa and 
the rotation number around main islands of pa is 1/9, 
that is not in the golden mean sequence. This orbit is lo¬ 
cated in the center of stable islands which are embedded 
in the chaotic layer around principal resonances. This 
orbit also does not exist at AT = 0. Generally, the anal¬ 
ysis performed on the two symmetry lines confirms the 
self-similar structure of diffusion and demonstrates the 
important role of non-standard resonances. 

The above analysis allows to understand some impor¬ 
tant properties of local diffusion rate D near the critical 
invariant curve Tg. The self-similar structure of D shows 
also the importance of non-standard resonances, different 
from principal approximants r„ = Pn/qn of Vg. These 
non-standard resonances are also self-similar at differ¬ 
ent renormalization levels. However, their sizes in phase 
space are comparable or sometimes even larger than those 
of principal resonances. Therefore trajectories can be 
trapped for a long time around these non-standard reso¬ 
nances and diffuse very slowly to internal chaos bound¬ 
aries surrounding islands of these resonances. Since the 
sizes of non-standard and principal resonances are com¬ 
parable, the contribution of internal chaos boundaries to 
Poincare recurrences may be relevant, and thus the de¬ 
cay of P{t) may be non universal. The determination 
of the asymptotic behavior of P{t) requires also a better 
understanding of transition rates between different renor¬ 
malization levels, which are not directly related to local 
diffusion and should be studied in more detail. 
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FIG. 1. Local diffusion rate D near principal resonance 
convergents r„ = Pn/qn to Vg = (5^^^ — l)/2. Do = K^/{8n^)-, 
chaotic component near unstable points with a) T = 10® 

V™ = 10, Np = 10 (x); b)T = 10^g„, V™ = 100, Np = 100 
(0)i cjregular component near stable points with T = 10®g„, 
Vu, = 10, Np = 10 (triangles). The straight line shows the 
theoretical slope a = 5. 


FIG. 2. Renormalized diffusion rate Dr/D o vs. rescaled 
position Aj/ij in the m—th level interval of renormalization 
scheme on the unstable symmetry line x = 0.5; m = A (full 
curve) and m = 6 (dashed line with crosses). The letters in¬ 
dicate the drops of diffusion Dr (see text). Other parameters 
are as in Fig. 1 a). 


FIG. 3. Renormalized diffusion rate Dr/D o vs. rescaled 
position /S.yR in the n—th level interval of renormalization 
scheme on the stable symmetry line x = 0.0; n = 8 (full 
curve) and n = IQ (dashed line with crosses). The letters 
indicate the drops of diffusion Dr (see text); T — 2 x 10^g„, 
V™ = 10, Np = 10. 
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